Searching for Secluded Dark Matter via Direct Detection of 
Recoiling Nuclei as well as Low Energy Electrons 
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Abstract 



Motivated by recent cosmic ray experimental results there has been a proposition for a 
scenario where a secluded dark matter particle annihilates, primarily, into Standard Model 
leptons through a low mass mediator particle. We consider several varieties of this scenario 
depending on the type of mixing among gauge bosons and we study the implications in novel 
direct dark matter experiments for detecting low energy recoiling electrons. We find significant 
CLY event rates and time modulation effects, especially in the case where the mediator is massless, 

that may be complementary to those from recoiling nuclei. 

1 Introduction 

> " 
00 . 

The analysis of the positrons excess (vs electrons) seen in cosmic ray spectra from PAMELA [1,2] 
in the energy region above 10 GeV confirming previous results from HEAT [3,4] and AMS-01 [5] 
experiments together with very recent results from FERMI [6] and HESS [7] collaborations seems to 
suggest the presence of a WIMP that annihilates into leptons without any indication of annihilation 
ON \ into (p,p) pairs or other hadrons (see Refs. [8,9] for relevant analysis). This is also reinforced by 
ATIC [10] experiment which reports excess of electron plus positron cosmic ray events in the energy 
region 300 < E < 800 GeV and also by signals from WMAP and EGRET [11-13] experiments. 
These phenomena can be explained by a scenario, originally proposed in ref. [14] - a subset of the 
so called secluded Dark Matter scenarios [15] - involving a new gauge boson xj}\, which couples 
to Standard Model (SM) particles and the WIMP through kinetic vector boson mixing with the 
following properties [17] : 

2m e <m x < m x (3 < m x a DM , (1.1) 

where m x /3 is a typical non-relativistic WIMP momentum and velocity j3 ~ 10~ 3 inside the galactic 
halo and «dm is the dark matter coupling. It has been shown that if eq. (jl.ip is satisfied then dark 
matter annihilation inside the halo to leptons is enhanced by a Sommerfeld factor of 0(«dm//5) [18] 
while annihilation to protons is simply kinematically forbidden. A typical range of parameters that 
are going to be exploited in our analysis and satisfy eq. (jl.ip are : mx = 0.1 — 1 GeV, m x = 0.1 — 1 



1 In earlier models [16] of secluded dark matter, WIMPs could be annihilated into new light scalar and gauge 
bosons. 
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TeV and odm = ct cm . The new force mediated by the X-boson is a long range force indeed. We 
must note here that there is a choice of another viable possibility with an even lighter mediator 
in MeV range that has been studied in ref. [19]. Our results for detecting low energy electrons are 
even more pronounced in this case. 

There is also a possibility for the gauge boson mediator to couple to the SM gauge bosons 
through a mass mixing matrix in a generalized gauge invariant way. These models are frequently 
called Stiickelberg models [20, 21] and are denoted as model type II in our classification. A char- 
acteristic of these models is that the electromagnetic current couples to the dark sector through a 
massless pole identified as the physical photon. As we shall see, this results in considerable and 
comparable rates in both nucleon or electron recoiling experiments. 

Alternatively, it could be that there is a symmetry that renders dark matter particles lepto- 
phylic [22-26] . This symmetry is spontaneously broken resulting in a massive gauge boson that 
couples directly to both leptons and WIMP at tree level. Again Sommerfeld enhancement dictates 
the mass of the X-boson to be in the GeV (or sub GeV) range. This is the model III that we 
consider in Chapter 2. 

Within the three model categories mentioned above we want : 

1. to study the implications of this new force carrier on both traditional nucleon recoil, and 
untraditional electron recoil direct dark matter searches, and, 

2. to suggest new dark matter experiments involving the detection of electrons scattered by 
this carrier providing a direct link to the recently observed cosmic ray anomalous elec- 
tron/positron events. 

So far there is a dedicated analysis for electron recoils in DAMA experiment [27] with energies 
approximately 5 KeV. Our analysis investigates recoiling electrons with energies as low as 10 
eV, and suggests an experimental method on how to reach such low energies. It is therefore 
complementary to the analysis of Ref. [27]. 

The structure of this article is as follows : In section 2 we present a field theory setup which 
helps to categorize three representative model examples that have recently been studied in detail. 
In section 3, we present event rate predictions for conventional nucleon recoil detection for the 
models studied. In section 4, we deal with the not so familiar methods of electron recoil detection 
rates together with time modulation effects. We also make a proposition of a prototype experiment 
to be exploited in discovering low energy electrons ejected from WIMP + atom collisions. In section 
5 we present our conclusions. 

2 Theory Setup and Model Categories 

In this section, we formulate the problem of the Standard Model coupled to, for simplicity, an 
Abelian dark sector with arbitrary kinetic or mass mixing terms allowed by Lorentz, gauge sym- 
metries and renormalizability. Our formulae are then applied in subsequent sections to make 
predictions for event rates in dark matter detection experiments. 

To read out the gauge boson propagators we start by writing the general renormalizable form 
of the Lagrangian : 

Jgf = -^S^ZC*"" + ^*J.M 2 ^ - ^d»&l~d v &„ + J/* M , (2.2) 
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where <& M j, = {d^ v — d u &^) is a A^-column matrix field strength tensor corresponding to a N- 
column <fr^ vector field, "T" denotes the transpose of a matrix, K and M 2 are real and symmetric 
N x N matrices with model dependent elements to be specified below and 3 is the gauge fixing 
N x N symmetric matrix necessary to remove unphysical gauge degrees of freedom. Interaction 
terms are encoded in the last term of eq. (|2.2p where an external current J M associated with 
symmetries, couples to the gauge fields. 

One has to notice that elements of the mass matrix M 2 should be further restricted by elec- 
tromagnetic gauge invariance. Phenomenologically speaking, there should always be a pole on 
the propagator ($>^& u ) corresponding to the massless photon i.e., the determinant of the inverse 
propagator at zero momentum must be exactly zero. Furthermore, without loss of generality, we 
can always assume that the diagonal elements of /C are normalized to unity. 

It is standard textbook exercise to find the Feynman propagator, T> fiu (p) with momentum p, 
for the gauge field 3?^ which in momentum space reads, 

iv^ip) = (iCp'-M 2 )- 1 (<v-^r) + W-M 2 )- 1 Ml. (2.3) 

At lowest order in K, interactions among fields are stored in the action functional 

S[J] = lj JJ(P) J*(-P) > ( 2 - 4 ) 

where J^(p) is the vector current in momentum space. Eqs. (|2.3p and (|2.4p are what we actually 
need to describe observables that arise from mixing dark (or hidden) and visible gauge bosons. As a 
simple example, consider the electromagnetic and the dark gauge boson current. Then in eq. (|2.2p . 
it is = (eJ^ ,m , gxJn '• It is then clear from eq. (|2.4p that interactions between the visible 
and the dark sector will involve off diagonal elements of the propagator (12. 3p . Observables, like 
nucleon recoil event rates can easily be described using the above propagator mixing formalism [28], 
by simply finding the inverse matrices such in eq. (|2.3p for a given model. We remark here that 
the propagator mixing formalism works equally well in different current basis such as Q — T3 or 
Y-T 3 . 



2.1 Model I : Non-standard Kinetic Mixing K, 

Models in this category [14,15] have been recently exploited in ref. [17] as candidates for explaining 
positron excess in cosmic ray data experiments. In its simplest form, the dark matter particle, 
X, is charged under a 'dark' U(l)x and the corresponding 'dark' gauge boson mixes with the 
photon Aa and Z-gauge boson, Z^. Annihilations of dark matter particles into only SM leptons 
(and not quarks) are kinematically allowed when the intermediate gauge boson has a mass at the 
GeV scale. 

In notation of ref. [29] and in basis (A^, X^, Z^) (or else Q — T3) our matrices K. and Ai 2 
appeared in eq. (|2.3p . become: 

\ / \ 

esin^H/ , M 2 = m\ , (2.5) 

1 / \ m| / 

where mx is the mass of the exotic gauge boson, mz is the mass of Z-boson, 9w is the weak 
mixing angle and e is a small (~ 10~ 3 ) mixing parameter between U(l)y and U(l)x field strength 



K 



1 

-ecos 8\y 




-e cos 6w 
1 

e sin 6w 
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V\AA/g)V\AA/ J = & , . (p) 

Figure 1: Diagramatic form of Feynman propagator appeared in eq. il2,3|) between gauge boson 
''flavours" i and j. For explicit expressions in model I see Eqs A2.6}) - i2.11\) : for model II see 

eq. mm . 



tensors. Working in Feynman gauge (H = 13x3) and keeping up to e 2 -terms it is easy to work out 
the mixed propagators V^ u (p), depicted in Fig(TJ between photon, X and Z-gauge bosons, labeled 
1,2,3, respectively : 

i5» = % + +0(e »), (2.6) 

■-fin, ^ ecos6 w f P^pA nr 3, , n „s 

p pZ _ y p J 



*2?^(P) = - ( o_^2 v^- rn 2\ \9^-^r) + 0(e), (2.1 



e 2 p 2 cos sin 6w ( P^Pv 
(p 2 — m 2 x )(p 2 — vn? z ) \^^ v p- 

VJ ,,.„ _ <W , E V (P 2 - cos 2 9„.mi) / „„,,„ 



i^tp) = + , 5 ^ 2 , 2 (a.,-^) + o( £ 3 ). (2.11) 

p pz _ m ^ i^pz _ rrr )yp — m z ) \ P J 

Some remarks are in order : i) among the three physical masses only m\ mass is shifted by 
an amount of m 2 x e 2 that we ignore ii) gauge invariance for the off diagonal propagator terms is 
preserved as should be the case. As far as the effective action eq. (|2.4|) is concerned, additional 
statements are in order: 

• The single pole [l/p 2 ] appears only in J e m • J e m exchange as usual in the SM. 

• A pole [l/(p 2 — m^)] for the exotic boson appears, apart from Jx • Jx exchange, also in 
Jem ■ Jx exchange at O(e). 

• There is exchange of current Jx • Jz i- e -> neutrinos and dark matter particles, through a 
double pole of X and Z at order e. 

• There is exchange of J em • Jz at order e 2 via double pole of X and Z 

The e ~ 10 _3 -term in the kinetic mixing can naturally arise as a result of mixing two U(l)'s at high 
energies - a mechanism that it was first proposed in Ref. [14]. Furthermore, X-boson contributions 
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to the muon anomalous magnetic moment relative to the SM expectation, Aa^ = a^ xp — = 
(290 ± 90) x 10~ n [30], are easily found using eq. (EHJ) to be 

Aq„ = ^ e 2 cos 2 9 W ( 2 , for ^ « 1 . (2.12) 
3ir \mx J nix 

This requires e<3x 10~ 2 for mx — 1 GeV where the equality accounts for the 2a upper limit on 
Aa^. Of course there are many other constraints on the mixing parameter e from direct or indirect 
collider searches and we refer the reader to recent work in Refs. [31-34]. For example, as we see 
from eqs. (I2.6p . (12. 9\} and (12. lip corrections to oblique electroweak observables arise at order e 2 
similar to the case of muon anomalous magnetic moment. 



2.2 Model II : Non-standard Mass Mixing, M 2 

Models belonging to this category are usually referred to as Stueckelberg models [20]. A recent 
account on "Stueckelberg" extensions of the Standard Model can be found in Ref. [35] . Here, it is 
more convenient to work on Y — Yx — T3 basis (B^, X^, Afy. We now assume that only the matrix 
Ai 2 is nontrivial, 




QyV 2 + my my mx 
2 



my mx 
-\ 9Y9 v 2 



m 



x 








(2.13) 



where gy,g are the U(l)y, SU(2)l gauge couplings respectively, m\ is a mass term for the hyper- 
charge gauge field and v is the vacuum expectation value. The form of the upper left 2x2 
Ai 2 matrix guarantees electromagnetic gauge invariance i.e., massless photon. Furthermore, the 
zero elements (23) and (32) guarantee that neutrinos are not charged under electromagnetism. 
Demanding that the inverse propagator has poles at the physical masses, det[p 2 — Ai 2 ]\ p 2 =m 2 = 
where = 0,mx,mz, we find that the photon mass is zero to all orders in my, the dark gauge 
boson and the Z-boson masses are not altered up to 0(m Y ), and thus m? z = \{g 2 + 9y)v 2 + 0(m Y ) ■ 
Following eq. (12. 4p we obtain the following effective action, 



S[J] 



d A p 
f2vr)4 



Je.m(p) ■ Je.m(-P) 



„ 2 9X my 

2e Je.m{P) ■ Jx{-P) 

gy m x 



1 

p 2 



+ 



9x Jxip) ■ Jx(-P) 1 



in 



X 



VTLr. 



+ 2e 29 -^^J e . m (p).J x (-p) 
g Y m x 



mxmy 

Zgvgx 5 — J x\p) • Jy{-p) 

m% 



p 2 — m 2 x 



+ 



m 



9 Jzip) ■ Jz(-P) + 9x^rJx(p) ■ Jx(-P) 



rrir 



+ igygx mx ™ Y Jxip) ■ Jy(-p) 



VTLr. 



2 1 2 \ + 0(m 2 Y ) , 
p z — m z 



(2.14) 



where e = gyg/\J 9y + g 2 is the electron charge. Furthermore, J e . m (p) = Ja 3 (p) + Jy(p) is the 
momentum space Fourier transform of the electromagnetic current, i.e., Je.m = Y2f Qffl^S with 
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Qfe being the charge of a generic fermion /. The dark current Jx obtains an analogous formula 
with obvious replacement of charge Qfe by another (hyper)charge, Qx- Of course, if fermions 
under consideration are Majorana particles then the corresponding current has only axial-vector 
form. In addition, Jz denotes the Fourier transform of the Standard Model neutral current = 
*<) (Ja — sin 2 6 W Je.m) where the electromagnetic current is, as usual in the SM, the sum of the 

COS u UJ x -/13 

third component of the isospin and hypercharge currents Jy. 

The physics of eq. (|2.14p is now transparent : to order O(my), there are interactions between 
the electromagnetic J e . m and dark current Jx mediated by the photon i.e., the dark matter particle 
is charged, and interactions between the hypercharge Jy and dark current Jx mediated by (X or 
Z) gauge bosons, respectively. An estimate of the dominant contribution to Aa^ results in an 
upper bound <9x 1CP 4 , where a 2a bound on Aa M is taken from Ref. [30]. 

2.3 Model III : Direct coupling, no mixing 

In this model, some of the SM leptons (but not quarks) £l,£r and the WIMP particle x are coupled 
directly to the dark gauge boson X^ in principle with different coupling^! : 

J x = g'Y\e L )l Ll n L + gfY'{e R ) en-fen + gxY'( x )x^X, (2-15) 

where Y'{ei,,en) = (lj — 1) denotes the particle hypercharge under the new gauge symmetry. As 
it has been suggested in Ref. [8,22,23,25,26], this could be an anomaly free gauged U(1)l,,~l t - 
Of course a new Dirac fermion x would be playing the role of dark matter particle is also gauged 
under this symmetry with Y\x) = 1- Because we have already discussed the effects of the kinetic 
and mass mixing in the previous models, without loss of generality, we assume that these mixing 
matrices are trivial in this model at tree leveH. If X^ does not couple to the muon then the most 
important constraint on a' = g' 2 /47r will arise from the v — e scattering at low q 2 : 

4" £ 7 x KT 7 . (2.16) 

m x 

We shall use this bound when discussing electron recoil detection rates in section 0] as is typically 
comparable (most of the time better) with other direct experimental bounds arising from LEP 
or meson factories. If the X^ vector boson couples to electrons and muons instead then there 
is a comparable bound to eq. (|2.16|) from the muon anomalous magnetic moment. Following 

Aa u = % — t- for m u <C my, there is a bound 

a' 

—5- < 4.4 x KT 6 . (2.17) 

m x 

3 Conventional WIMP searches 

Conventional DM searches deal with phenomena of WIMPs scattered of a nucleus. The study of 
the recoil energy spectrum is the primary goal of experiments such as CDMS [36], XENON [37] 
and DAMA [38]. For models we described in the previous section there are two cases which have 
been discussed recently in the literature that could explain the recent anomalous cosmic ray events: 

2 Various possibilities on how this is realized can be found in Ref. [22]. 

3 Of course mixing of the gauge boson with the U(1)y is inevitable at one loop. Its magnitude is calculable : 
e~a 2 log^=2x 10~ 4 for a' = a cm . All the rest will then proceed following eqs. (|2. 61 - 12. lip of model I. 
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Figure 2: A Feynman diagram leading to the direct interaction of the WIMP \ to the quarks 
relevant for direct detection of dark matter. The process is mediated by the physical photon. The 
cross indicates merely that the exotic gauge boson has a small admixture of the photon. Similarly 
the WIMP can also couple to electrons. 



a) The lightest mediator is massless and 

b) the lightest mediator is massive with mass around the proton mass (m p ), 
in addition to the assumption that 

m p <C m x , (3.18) 

where m x is the WIMP mass. Only model II belongs to the first category and models I, II belong 
to the second since by definition, there is no direct coupling of X-boson to quarks in model III [see 
however footnote 2]. In the following subsections we present the WIMP-nucleon cross section for 
both cases (a) and (b). 



3.1 Massless Mediator 

The differential WIMP-proton cross section in the rest frame of the initial proton is given by: 

s(f3) e 2 (gxn) 2 d 3 p' d 3 g 3 (3) , 
p -a J^rf (2^)3 ( ] 6 (P " P " q) ( } 5{ ~ ~ q) ■ (3 9) 

In the above equation p', p are the momenta of the initial WIMP and the final WIMP and q 
the momentum transfer to the nucleon and T = p 2 /2m x , T' = (p 1 ) 2 /2m x and T q = q 2 /2m p , are 
respectively the corresponding kinetic energies in the non relativistic limit. Furthermore, (3 is the 
WIMP velocity and s((3) = 1 for a WIMP which is a Dirac fermion, while s((3) = 1 in case it is 
Major ana one [39] One finds that the momentum transfer and the final nucleon energy are given 
by: 

q = 2p r v^ « 2m p v^ , T q m 2m p v 2 ^ 2 , (3.20) 

where \i r is the WIMP-nucleon reduced mass, m p is the proton mass and < £ < 1 is the cosine of 
the angle between the incoming WIMP and the outgoing nucleon. Integrating over the momentum 
of the outgoing WIMP and the magnitude of the momentum of the final hadron as well as the 
0-angle one finds : 

^fM'*. (3. 21) 

2vr (2m p ) 2 v 3 £ 3 v ; 



4 The Majorana fermion does not possess electromagnetic properties. Hence only the 7 M 7s of the WIMP -X-boson 
interaction contributes. 
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The above expression exhibits, of course, the infrared divergence. We will impose a low momentum 
cut off E^/ A provided by the energy threshold where A is the mass number of the target, i.e. 



Un ~ \/ (2Am p f3 2 ) ■ (3 - 22) 
Thus the total cross-section for a Majorana WIMP is given by: 

a , 2 1 f Am p m p \ a , 2 1 Am p 



a = - (g X KY — ^ — ^ - « - ^2 ^ • ( 3 - 23 ) 

2 (m p ) 2 V ^th Jmax/ 2 (»™p) .E th 

Eq. ()3.23p shows a much stronger dependence of the event rate on the threshold energy i? t h due to 
the adopted cut-off -^cut-off = Eth/ A. It is interesting to note that this cross section is independent 
of the WIMP velocity (in the case of a Dirac WIMP the extracted from the data cross section must 
be multiplied by (3 2 ). We distinguish two cases : 

1. The case of Majorana WIMP. We find: 

a « 1.6 X 1(T 30 cm 2 (g xK ) 2 . (3.24) 

Eth 

The direct dark matter experiments have recently set on the coherent nucleon cross section 
the limits: 

• The CDMSII experiment [36]: 

The best limit is 6.6 x 10 -44 cm 2 . The extracted value depends, however, on the assumed 
WIMP mass. So it can vary between 6.6 x 10 -44 and 6.6 x 10 -42 cm 2 . 

• The XENON10 collaboration [37] 

They extract 8.8 x 1CT 44 cm 2 and 4.5 x 1CT 44 cm 2 for WIMP masses of 100 and 30 GeV 
respectively. 

For our purposes we will assume that the extracted from the data nucleon cross section 
is 10 _7 pb = 10 _43 cm 2 . Furthermore we will take as a reference a threshold energy of 5.0 
KeV and examine the sensitivity of our results to the experimental threshold. Using the 
experimental limit, a v < 1.0 x 10~ 43 cm 2 , we can write: 

Rate (new) R Z 2 . , 9 Am„ 

-, i 1 — - = 1.6x10^ 3xk 2 —^. (3.25) 

Rate(conventional) A z Eth 

Note that the coherence factor now is Z 2 , since in the case of the photon only the protons of 
the target contribute. Adopting a threshold value of 5 KeV, we get 

Rate(new) ^ 18 Z 2 2 (3 26) 

Rate(conventional) A 

For the Ge target (A = 73, Z = 32) we get 

4.3 x 10 iy (gx^f , (3.27) 



Rate(new) 



Rate(conventional) 



8 



which leads to the limit: 



\9XK\ < \ 77T = 1-6 x 1(T 10 . (3.28) 

|y 1 ~ V 0.43 x 10 19 v ; 

Prom the second term in eq. (I2.14p and assuming that ckdm = 9x1^ = a em one can easily 
translate this into bounds on the model II parameters for Majorana WIMP : 

Qx — < 0.54 x 1(T 10 , model II . (3.29) 

m x 

2. The case of a Dirac WIMP. We find: 

~ 1 a 1 ( \ 2 Am P a qpa 

■ (3 - 30) 

If we knew the coupling \gx&\ we could incorporate this into the evaluation of the nuclear 
cross section, fold it with the velocity distribution and proceed with the evaluation of the 
event rate. Since, however, we like to constrain the parameter \gx*\ we will employ an 
average velocity: 

°^<°>™<j.> 2j^r {9xK) • (3 - 31) 

But for a Maxwell - Boltzmann distribution i.e., < > —> < ^ > , we obtain the constraint: 



\g X K\ < 1.6 x iQ-io V^<^ 2 _^ ra 0>8 x 10 -i3 ^ (3.32) 



from which the bound on model II for «dm = a em , 

Qx — < 0.27 x 1(T 13 , model II , (3.33) 
is found. As expected the limit is now more stringent than in eq. (|3.29ft . 

The results for the Xe target are similar. This bound is by many orders of magnitude stronger 
than the one obtained from electroweak fits [35] or {g — 2)„ [see discussion towards the end of 
section (|2.2j) ] . The corresponding bound for Dirac WIMP is about three orders of magnitude more 
stringent. This means that additional mechanisms should be added in model II (Stuckelberg type 
of Ref. [35] for example) in order to efficiently depleting the WIMP in the early universe (the 
diagram in Fig. [2] is just the crossing diagram of the annihilation cross section). 

Although eq. (|3,29[) [or eq. (|3.33|) ] provides a very stringent limit, we should not forget that in 
this case we have a much stronger dependence of the rates on the energy threshold through the 
need for a low energy cut off on the elementary cross section. 

Alternatively we may extract from the data for Xe (A=131,Z=54) an elementary cross section 
assuming it to be of the fornH : 

<? / \ \ -4 5 keV 

a s N>x0 (A,E th ) = a — , (3.34) 



5 This treatment does not distinguish between a Majorana and a Dirac WIMP. 
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Figure 3: The total rates for traditional WIMP searches assuming a nucleon cross section cj at = 
10 3 cm 2 in (a). The case of the photon mediated process considered in this work is exhibited 
in (b). Both refer to the case of a heavy target (A=131) and were computed assuming an energy 
threshold of 5 KeV. The results for the Iodine target used by the DAMA experiment are almost 
identical. 
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(a) (b) 
Figure 4: The same as in Fig. \3[in the case of the light target 19 F. 



where gq is the elementary cross section obtained in the particle model for a target with nuclear 
mass number A and threshold energy E^. Then by fitting to the experiment we obtain 

(131/54) 2 ct^ x o = 0.5 x 10" 7 => ct = 2.9 x 10" 7 pb = 2.9 x 10" 43 cm 2 . (3.35) 

In spite of the (Z/A) 2 factor we obtain a smaller value than in the standard experiment. This is 
due to the small cut off energy E^/A employed. With the above ingredients the number of events 
in time T due to the coherent scattering [40], can be cast in the form: 



-5 

" ' '' ''N, x 



T p(Q) 100 GeV m of* 

ly 0.2 GeVcm -3 m x o Tkg 280 km s _1 10" 43 cm 2 



f coh (A,fi r (A)) , (3.36) 



where the elementary cross section afj x can be treated as a phenomenological parameter indepen- 
dent of the WIMP mass in units of 10 -43 cm 2 . The quantity f co h(A,/j, r (A)) can be obtained from 



10 




Figure 5: The quantity R(E t ^) / 'R((E t h) min , i.e. the ratio of the event rate at a given threshold 
divided by that at the lowest threshold considered, as a function of the threshold energy. In (a) as 
predicted by traditional mechanisms (lowest threshold assumed zero). In (b) as predicted by the 
present model (now due to the need for a cut off the lowest threshold energy employed was 5 keV). 
The thick line, short dash, long dash, hne line and long short dash correspond to WIMP masses 
10, 50, 100, 200 and 500 GeV respectively. 



the published in Ref. [40] values of t for the standard MB velocity distribution (n=l). For the 
photon mediated mechanism examined here the above equation must be modified by multiplying 
fcoh{A, /x r (j4)) with the factor Z 2 /A 2 and employing eq. (13.340 for the elementary cross section (in 
units of 10 -43 cm 2 ). The event rate per kg of target per year for the traditional experiments for a 
heavy isotope like Xe and a light isotope like 19 F, as a function of the WIMP mass is exhibited in 
Figs [3] and 01 On the same plots we show the event rate for the photon mediated process examined 
in the present work. It is not surprising that the agreement is good since the elementary cross 
section was fitted to the data. The small difference is understood, since in the extraction of the 
elementary cross section from the data a zero threshold value was used in the phase space integrals. 
The event rates are sensitive functions of the threshold energy, R = R(E t h). In the case of the 
Xe isotope the ratio R(E t h) / R((E t Y 1 ) min is exhibited in Fig. [5j The threshold dependence is much 
more profound in the case of the light WIMP, since, then, the average energy transfered is small. 
As expected the threshold dependence is more dramatic in the case of the present model (this is a 
bit obscured in the figure since in this case the graphs are normalized at 5 keV). 

In the case of a Dirac fermion the extracted limit will be smaller, but the traditional calculations 
are not adequate for the analysis, due to the different velocity dependence of the elementary cross 
section. 



3.2 Massive Mediator 

In this case the WIMP - nucleon cross section reads : 

167ra em k 2 aj)M m 2 
a = s{f3) 

m x 

_ u^^^ ( , 37) 
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where the cross section refers to Dirac (Majorana) WIMP and s(/3) = l(/3 2 ) respectively. Taking 
P 2 -> < P 2 > w 10" 3 we find: 

k<3x 10~ 7 (3 x 10~ 4 ) . (3.38) 
From these we obtain bounds for parameters in models I, II [see eqs. (I2.8P and (12.140 ] , 

e < 3.0 x 1(T 7 (3.0 x 10~ 3 ) , Model I (3.39) 

Qx—< 1.6 x 10~ 6 (1.6 x 10~ 3 ) , Model II , (3.40) 

m x 

where the number in parenthesis corresponds to Majorana WIMP dark matter particle. These 
limits are less stringent than those obtained in the case of the massless mediator. 
In the case of the massive mediator, with the possible exception of the velocity dependence in the 
case of Majorana WIMP, the cross section behaves as in the standard CDM case, since in this case 
we do not encounter an energy cutoff. Since, however, we do not know the values of the parameters 
e and we cannot make predictions about the event rates. Instead we have used the present 
experimental limits to constrain these parameters. Thus we saw that the current experimental 
limits impose the most stringent limits on these parameters. If, on the other hand, we use the 
previous constrains we can conclude that WIMPs in models I, II scatter off nuclei too many times. 
These effects should have been seen in experiments [36,37] (or may have already been seen [38]). 
An exception is a Majorana WIMP candidate in model I which results in current sensitivity event 
rates. 



4 Unconventional WIMP searches 
4.1 Cross Section 

The other possibility is the direct scattering of WIMPs by electrons that are bound in atoms. The 
relevant Feynman diagram is depicted in Fig. [2] with quarks replaced by electrons. In this case 
only the electron flavor can be detected since the other flavors are not energetically allowed. Since 
the outgoing electrons are expected to have energies in the eV region one cannot ignore atomic 
binding effects. The binding energy b is found from the tables of ionization potential (energy) of 
an atomic 

The problem is to find the cross section for WIMP scattered off an electron bounded in an 
atom. In order to proceed we shall make two simplifying assumptions : 

1. As a working example, we shall assume that the target is a hydrogenic atom denoted by H 
i.e., a nucleus with charge +Ze and a single bounded electron with charge — e. We shall 
discuss deviations from this assumption throughout. 

2. The gauge boson mediator X couples only to WIMP and leptons but not to quarks. This is a 
necessary condition to explain PAMELA positron excess of events. Therefore, this discussion 
refers strictly to model III in eq. (|2.15|) [see however footnote 2]. 

6 Tables are normally given in kj/mol, but they can easily be translated in eV, since 96.485 kj/mol = 1 eV. Thus 
for Cs we find b = 375.7/96.485 = 3.89 eV. 
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There are four processes that could take place in WIMP + H-like atom collisions 



X + H — ► x + H (elastic) , (4.41) 
X + H — ► x + H* (inelastic) , (4.42) 
X + H — > x + e~ + H + (production) . (4.43) 

For the rest we shall consider only the situation (|4.43|) . The elastic scattering (|4.41|) cannot be 
detected, and although we cannot exclude the inelastic one (|4.42jl from being experimentally probed 
through final state photons, we believe that it would be easier to detect the electrons from (|4.43|) . 
We shall assume that the electron emerges with high momenta, p' e , such that in the final state 
its interaction with the Coulomb potential in H-like atom is negligible, i.e, we can use plane wave 
states for incoming and outgoing particles. Using standard textbook [41] wavepacket analysis our 
starting point will be the cross section formula in the lab frame: 

x d 3 p.(2^) 3 i< 3 )(p x + p,-p^-pa|(6(Z,p,)| 2 , (4.44) 

where p x , p e (p^, p' e ) are the incoming (outgoing) three vector momenta of the WIMP and electron 
particles respectively, and M is the matrix element of the process x + e ~> X + e averaged over the 
spins of the initial states calculated in Born approximation. We also ignore local velocity effects 
from the bound electron in the (static in lab frame) atom i.e., that is the relative velocity is v — v x . 
T{ = pf/2rrii,i = Xi e are the kinetic energies and b is the binding energy of the electron in H-atom 
(~ 13.6 eV). Moreover, in non-relativistic limit E x ~ E' ~ m x and E e ~ E' e ~ m e with m x 3> m e , 
while 4>ntm t {v)i normalized at J v d 3 p\4> n £ me (p)\ 2 = 1, is the Fourier transform of the coordinate 
wave function ijj n £ mi (r). Using the 5^ -function to perform the integration over p e , we obtain: 



. \M\ 2 d'Px^Pe JlPxl 2 IPxl 2 IPel 2 h(7 ,\ u f7n r M , 



)| 2 ,(4.45) 



where the energy conservation delta- function has been written out explicitly. The result of eq. (|4.45p 
is a product of two parts : a part that contains the dynamics of the WIMP-electron interaction 
through the matrix element \M\ times the probability of finding the target electron with momen- 
tum p e = p' x + p' e — Px m H-atom. In addition the matrix element of the process x + e — > X + e 
averaged over the spins of the initial states in Born approximation reads : 



(16-7T) 2 'a-EtMOt' 'rrilrri^ 
(IPX - Pyl 2 - m \f 



where the factor s((3) = 1 (/3 2 ) for Dirac WIMP (Majorana WIMP) particle. Note that the cross 
section for Majorana WIMP is always smaller by a factor of /3~ 2 compared to the one involving 
Dirac WIMP. We now use the kinetic energy 5-function appearing in eq. (14.450 in order to perform 
the lp' I integration and arrive at: 

d(J = S ^)l\ 7T2 TY2 r^lPel^lPel \<j>nim e (Z , p' x + p' e - p x )\ 2 d£ d V , (4.47) 

vlPx P%l m x> IPxl 
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where the initial WIMP momentum is |p-J = m x [3 and the scattering angles are defined as 

C = P x -p'x' ri = P x -p'e, Z,V G [-1,1] ■ (4-48) 

The integration over the azimuthal angles has been carried out trivially in eq. (|4.47|) and the 
momentum |p^| of the scattered WIMP is found to be 



|p'J = Jm 2 f3 2 - 2m x b(Z) - ^ p' 2 , with p' e = ^2m e E' e , (4.49) 

where b(Z) is the ground state energy for hydrogenic atoms is 

Z 2 e 2 Z 2 1 

b(Z) = — — = —m e a 2 m , a~ , (4.50) 

2a 47T 2 m e a em 

2 

in the approximation fi ~ m e where /i is the reduced mass, with. Q em — ^— « 1/137, m e ~ 0.5 
MeV and a = a>o ~ 0.5 A being the Bohr radius for Z = 1. Throughout this chapter, we are going 
to use the ground state momentum distribution of hydrogenic atoms which reads: 

2 3 / 2 (Za) 5 / 2 

(j) 100 (Z,p) = 2^2- ( 4 - 51 ) 

7ra (Z z + p z a z ) z 

Notice that since 0ioo(p) depends on \p\ 2 and therefore from the scattering angles r\ and £ and 
electron energy £g. A term in eq. (|4.47p . j^j = j^j, arises from the fact that we treated the 
H-atom as a brick wall potential. Had we not done so, the influence of the Coulomb potential on 
the emerging electron would not have been uniquely correlated to p^,p x and the back reaction of 
the proton should have been taken into account. 

Exactly the same result as in eq. (I4.47|) can be found by using simpler time-dependent pertur- 
bation theory for transitions to continuum in non-relativistic quantum mechanics [42]. In a more 
refined analysis however, when the recoiling energy is in the neighborhood of the binding energy 
of the atom one should take into account effects from the continuum hydrogenic wave functions 
instead of treating the final electron as plane wave. This analysis, though more accurate, is far 
more complicated and does not change the qualitative features of our results. It can be addressed 
in the future (together with other effects) especially if these kind of experiments become operative 
[see subsection 4.4 below]. 

We analyze below the corresponding cross sections for a massless and a massive mediator as 
we did in section [3] for the nucleons. 

4.1.1 Event Detection Rates 

In general for an atom, due to binding energy effects only the loosely bound electrons can contribute 
to the process (|4.43[) . So we will convolute the elementary cross section with the WIMP velocity 
distribution, which, with respect to the galactic center, we will take to be Maxwell-Boltzmann 
form: 

3/2 1 „„2 



Transforming this into the local coordinate system: 



P^(3(3 + p z = (3l3 + \j 2 < f 2 > z , p 2 ^ p 2 + 1 < f3 2 > +2(3cos(6)J^ < p 2 > , (4.53) 
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where 9 is the angle between /3 and z and (3g = y 2< g > is the sun's velocity with respect to the 
center of the galaxy and < (3 2 10 -6 . Then we obtain the local distribution of speeds fe(f3) 
relative to the detector to be 

3 V' 2 i -(rife + 2,3 CO! «"^5F +1 ) 



The integration over the angles of the distribution can be done analytically. In evaluating the 
rate one has to incorporate the oncoming flux. So, adopting appropriate normalization, in the 
convolution we introduce the factor l/y< /3 2 >• This way we find the rate to be proportional to : 

(3MP)d(3 = ( 3 \ 3 / 2 2_ -(stfe+i) P 3 sinh (Wgg < /3 2 >)) 



_ e v 2< ' 3 > 7 — v — dd (4 55) 

y^T^ V2</3 2 >7 y<^^ /V 3 /( 2 < /? 2 >) 

Combining this with the cross section of eq. (|4.47j) obtained previously we arrive at: 

/ da fS \ ft>~ DS,{f>) da 

where the lower velocity in the integral can be read from the positivity of the square root quantity 
in eq. (gggj) 

_ /as + »w i (4 . 57) 

and /3 eS c = 2.84-^/ (2/3) < j3 2 > is the escape velocity. It is now easy to calculate the differential 
event rate per eV ejected electron energy per year and per kilogram of target material to be 

where p Q = 0.2 GeV/cm 3 is the WIMP energy density and N e is the number of target electrons. 
Integration of eq. (14.580 upon E' e over the region from E' e . = to [m x /3 2 sc /2 — b(Z)] results in the 
total event number per unit time and mass of the target which among other parameters depends 
on the mass and atomic numbers of the target atom. Moreover we shall display results on the total 
event rate R{Z) when E' e . = with varying experimental threshold energy E t ^. 

4.1.2 Time Modulation Effects for Electrons 

In the convolution of the elementary cross section we have so far considered only the motion of the 
sun with respect to the center of galaxy. More realistically, one should consider also the Earth's 
velocity and then find the modulated event rate that might be detected on Earth. In this case the 
WIMP velocity is read from 

v' = v + vq z + fi(sinax + cos a cos 7 y + cos a sin 7 z) , (4.59) 

where t>o is Sun's velocity, v% is Earth's annual velocity, 7 = ^ is the angle between the projection 
of vector vi on the plane yOz and the y direction and a = a(t) is the complementary angle of the 
angle between vi and x. Then the WIMP cross section has to be convoluted with 

'Pf e (J3)dp\ f[3Mf])df3\ , 

' — ' 1 (1 + fcocosa), (4.60) 
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where the expression with the subscript "0" refers to eq. (|4.55|) with 8 = ^ ~ 0.135 and 



cosh 



2 < > sinh (2/3^=1^ 

It is now trivial to extend the distribution with energies event rate of eq. (|4.58p with 

dR I dR\ / dR\ fl . 

J) + (l^) X«*« (4.62) 



dE' e \dE'J \dE e/mod 

where (^j§r\ is the unmodulated differential event rate while (-^r\ contains also the factor 
k in eq. (jlTBTj) . 

4.2 Massless Mediator 

In this case dark matter scattering happens via the coupling of the exotic gauge boson to the photon 
(model II). The Feynman diagram is identical to the one presented in Fig. [2] with the quark replaced 
by electron. In general case the WIMP-electron cross section is not independent of the velocity. 
Thus, we will first estimate the cross section by using an average velocity \/ '< f3 2 > = 10~ 3 . 
Following eq. (|4.47p for a photonic mediator we find the differential cross section: 

d ° . IO S^_%.,. ._2_2_ IPxl ._/ I f 1 ^ f 1 drj \<j>nem e (Z,P' X +Pe-Px)\ 2 ? 



IP' I f 1 f 1 

s((3)167r 2 a'a.T)M k 2 m 2 m e - — %j |p^| d£ 

\Px\ J -I J-i 



dE> v * C \P X \'^ J-i ^ J-i ' (P' X -P X Y 

where q = p^. — p x is the WIMP momentum transfer which is £ dependent. The cross section 
peaks up the most from the forward direction £ ~ 1. It should be mentioned that since the initial 
electron is bound, there is no infrared divergence in this case. Moreover, the momentum transfer 
can be as low as : 

ld.2M±S. (4 . 64) 

This relation is important for explaining our numerical results below. Furthermore, in presenting 
the results we assume a Dirac WIMP fermion,i.e. s(/3) = 1. Furthermore, we choose a benchmark 
scenario inspired by our findings in nucleon decay : 



(3 = v / < j3 2 > = 1(T 3 , Z = 1 , a DM = a' = a om , m x = 100 GeV , k = 10~ 10 . (4.65) 

As it is obvious from eq. (|4.63p it is very easy to apply our numerical results to any other parameters, 
/3, aDM) a', k, than those shown in eq. (|4.65|) . We must note here that there is no parameter k in 
model III. This parameter is used here as a rescale factor and its very small value is adjusted so 
that we obtain rates of few events. 

In Fig. are shown the results for the da/dE' e as a function of final electron's energy E' e 
for three different cases of hydrogenic atoms with Z = 1 , Z = 3 and Z = 6 respectively. The 
differential cross section takes on its maximum values for final electron energy of around few eV 
for Z = 1, around few tens of eV for Z = 3 and around a hundred eV for Z = 6. For the case 
Z = 1, the extremum happens because of a fast increase of the term j^j |p^| ~ \f~E'e an d the almost 
constant value of | <^>ioo 1 2 until 5 eV . For higher electron energies, e.g., E' e > 10 eV, the probability 
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Figure 6: a) Predictions for da/dE' e as a function of the ejected electron energy E' e . The target is 
assumed to be a hydrogenic atom in the ground state with Z=l,3,6 (from top to bottom), b) The 
total cross section for process &4.43\) as a function of the experimental threshold energy for two 
binding energies, c) The differential event rate as a function of the electron energy and various 
WIMP masses (10,100,1000) GeV from (top to bottom). Other parameters not shown, are taken 
from eq. (I4.6'5|) . 



density factor |0ioo| 2 drops fast as 1/E'^ 1 and the term in the denominator of the integral increases 
as .Eg 2 , resulting in overall decreasing of the cross section as E e 19 ^ 2 .The same analysis can be used 
to describe the behavior of dcr/dE' e in the other cases (Z = 3, Z = 6). We must note here the in 
the limit E' e — > we obtain da/dE' e — > as the case should be. This is obscured in Fig. [6] due to 
the range choice of E' e . 

Corresponding to the input parameters noted in (I4.65P we calculate the total cross section from 
eq. (|4,63p after numerical integration over E' e in the region [E t ^,m x (3 2 /2 — b(Z)]. Our results for a 
vs. the threshold energy £^h are depicted in Fig. [6b- We have chosen two extreme cases of binding 
energies : b = 0.74 eV that is the the binding energy of the electron bounded in the two electron 
atom H~ and b = 13.6 eV that is the one corresponding to the H- atom we have been dealing so 
far. For E^ < 10 eV the difference in cross section is about three to six orders of magnitudes while 
for higher threshold energies becomes unimportant. 

Following eq. (14.630 it turns out that the total cross section for process (|4.43j> is WIMP mass 
independent. It is experimentally useful to know how the cross section depends on the threshold 
energy E t h that a given experiment can accomplish. This is plotted in Fig. [6b- For E^ < 1 eV, 
the cross section is essentially independent of E t h- When the threshold becomes 5 eV, in the case 
of b(Z) = 13.6 eV, the cross section drops by a factor of 5 eV while up to 10 eV by a factor of 
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Table 1: Time modulation effects in case of a photonic mediator following eq. &4.62\) in the text. 
Various input parameters are given in eq. <I4.6'5|) . H is the ratio of the modulated divided by the 
unmodulated differential rate. 



50. For smaller binding energy though, i.e., b(Z) = 0.74 eV, and up to 10 eV the cross section 
decreases by three orders of magnitude. 

Furthermore, the dependence of differential event rate dR/dE' e as a function of the ejected 
electron energy E' e for three different WIMP masses, m x = 10, 100, 1000 GeV, is shown in Fig. [6b. 
There is a maximum which follows the behavior of differential cross section. The event rate falls 
as ^-/m x as the WIMP mass increases in accordance with eq. (I4.58p . For energy of few eV's and 
m x = 10 GeV we obtain a handful of events for k = 10~ 10 . A total event rate is obtained after 
integrating over the differential rate in Fig. [6b- As a typical value, for m x = 100 GeV and the 
parameters in (|4.65p we find R(Z = 1, k = 10~ 10 ) ~ 1 events/yr /target kgr. The reader must 
recall here that this assumes a mixing parameter as small as k = 10~ 10 !! 

Finally, following the theoretical discussion of the previous subsection we examine effects of the 
WIMP time modulation. In Table 1 we display both the unmodulated and modulated differential 
event rate for four representative values of E' e in the case of a massless mediator and parameters 
of eq. (|4.65p . The dimensionless parameter H, which is the ratio of the modulated by the non 
modulated differential amplitude, is constant around 9 — 13% independent of the energy and the 
WIMP mass. So the modulation h = 5 ■ k of the total rate is also going to be around 10%, which 
means that the difference between the maximum (here always in June 3rd) and the minimum (here 
always in December) is 18 — 26%, a result should not to be overlooked. 



4.3 Massive mediator 

Again the relevant diagram is the one of Fig. [2] with quarks replaced by electrons. By taking the 
non-relativistic limit of eq. (|4.47p and the assumption that the momentum transfer in eq. (|4.64p is 
much less than the mediator mass, q 2 ^ m|, we arrive at 

= s{(5) ^^m 2 m e -^-\p' e \ / di / dJ7|^ m ,(Z,p x + v' e -V x )\ 2 , (4.66) 

atj e m X IPxl J-l J-l 

In what follows we assume a Dirac WIMP fermion, i.e., s(/3) = 1. We assume the following input 
parameters : 

f3 = V '< /5 2 > = 10~ 3 , Z = 1 , a D M = ol = a em , 

m x = 1 GeV , m x = 100 GeV , k = 1 . (4.67) 

Although this parameter space violates the bounds in eqs. (|2.16p and (I2.17P it serves as a benchmark 
for this article in comparing results with those of section 3 if possible. The value of k is chosen 
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Figure 7: a) Predictions for da/dE' e as a function of the ejected electron energy E' e . The target is 
assumed hydrogenic atom with Z=l,3,6 (from top to bottom) in the ground state, b) The total 
cross section as a function of threshold energy, c) The total cross section as a function of mx for 
two different binding energies. We assume a Dirac WIMP, = eV and input parameters from 
eq. A4.67]) if not stated otherwise. 



such that the resulting rate presented in the figures assumes no mixing of the X-boson mediator 
which is formally the case of model III. 

Results for the differential cross section dcr/dE' e for the electron in the ground state of three 
hydrogenic atoms are shown in Fig. [Tfci. The differential cross section takes on its maximum values 
for final electron energy of around few eV for Z = 1, ten of eV for Z = 3 and around hundred 
eV for Z = 6. For the case Z = 1, the extremum happens because of a fast increase of the 
term j^j |p^| ~ s/E^ an d the almost constant value of |</>ioo| 2 until 5 eV [see eq. (|4.66|) ]. For 
higher electron energies, e.g., E' e > 10 eV, the probability density factor | <^>ioo 1 2 drops fast as 1/E' e 8 

' — 15/2 

resulting in overall decreasing of the cross section as E e . In physical terms, the outgoing 
electrons of high energy demand high momenta in the initial electron wavefunction, which leads to 
suppression. The dependence on the Z is easily explained if we recall that for hydrogenic atoms, 
(p 2 ) n= i = Z 2 Pq where po is the Bohr momentum for Hydrogen. Furthermore, despite appearances 
in eq. (|4.66p . the differential cross section depends only very mildly on the WIMP mass. One 
can show analytically that the double integral over the wave function squared, is approximately 
proportional to l/"i^ which cancels the m 2 , in the numerator. 

Corresponding to the input parameters noted in (I4.67P we calculate the total cross section from 
eq. (|4.66|) after numerical integration over E' e in the region [E th , m x (3 2 /2 - b(Z)] . For fixed velocity, 
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Figure 8: a) Differential event rate of Dirac WIMP scattered off hydrogen (Z=l, A=l) target 
electrons per year per Kgr as a function of ejected electron energy E' e in eV. Three different 
WIMP masses have been assumed : m x = 10, 100, 1000 GeV, from top to bottom, respectively, 
b) The total event rate as a function of the experimental threshold energy for m x = 100 GeV for 
two different binding energies. Other input parameters are taken from eq. i4.67\) for the massive 
mediator. 
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The total cross section increases by a factor of about 32 when [3 = (5 esc is taken. The cross section 
decreases with Z [see also FigJTK], the reason being the fact that the binding energy increases with 
Z 2 [see eq. (14.500 ] and therefore we need to go to larger - compared to ground state - momenta where 
the wavefunction is small despite their maximum value displacement towards larger momenta. 

Assuming that the sensitivity of detecting low energy electrons will be analogous to the ongoing 
experiments (~ 10 -43 cm 2 ), we could even extract bounds on various parameters in models I, II or 
III. From all running experiments, DAMA [27,38] is the one that triggers on final state electrons 
with energy around 5 KeV. From Fig. [7|i one obtains that, around that energy, the cross section 
is too small for mx = 1 GeV and all other inputs in eq. (|4.6Tj) . However, da/dE' e oc m x 4 and 
therefore for mx ~ 1 MeV i.e., model types proposed in ref. [19], DAMA is a relevant experiment. 
Additionally, this is demonstrated in Fig. [7fc where the total cross section as a function of mx is 
plotted for two reference values of binding energies. 

In Fig. [TJa we examine the total cross section as a function of the experimental energy threshold 
for low energies, relevant to our proposal. As we can see, the total cross section reduces by a factor 
of six in the region < Eth < 10 eV. Above 10 eV the cross section drops drastically [see total 
rate in Fig. [8b] . 

Although not shown, we have also examined departures of the wavefunction from the ground 
state. The maximum value da/dE' e \ max appears at the same place in E' e ~ 1 — 10 eV. As an 
example, the difference in da / dE' e \ max is an enhancement by a factor 20 when going from Is — ► 2s. 
Furthermore, the size of the momentum transfer in conjunction with the non-zero binding energy 
are such that never let the wavefunctions to reach their zero nodes. 
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Table 2: Time modulation effects in case of a massive mediator following eq. i4.62\) and various 
input parameters in eq. i4.67\) . H is the ratio of the modulated by the unmodulated differential 
amplitude. 



Assuming one electron per target atom, and the average cross section of Fig. [7h for Z = 1, 
the differential event rate per eV of electrons energy per year per Kgr of hydrogen material as a 
functions of E' e for various WIMP masses is depicted in Fig. [5k- The differential event rate again 
exhibits a maximum which follows that of the differential cross section calculated in Fig. [7Ji. The 
event rate is of course higher for smaller WIMP mass [recall eq. (|4,58[) ] and for electron energy of 
few eV's it varies from 0.01 up to 2 events/yr /kgr/eV for m x = 1000,10 GeV respectively. For 
electron energy of around 100 eV the role of the wave function is to reduce the differential rate by 
an order of magnitude i.e., from 10 -4 -j- 10 -3 events/yr/kgr/eV. The total event rate for m x = 100 
GeV and the other parameters in eq. (|4.67p is predicted to be: 

R(Z = 1,k = 1) ~ 2 [events/yr /target kgr] . (4.68) 

It is useful to know how the total rate (|4.68p varies with an experimental threshold energy. This 
information can be extracted from Fig. [5b for two different but judiciously chosen, values of binding 
energies. As in the case of the total cross section in Fig. [7b, the total rate drops by only a factor 
of five until E t ^ ~ 10 eV while it drops very rapidly after about this scale. For example, it drops 
by a factor of 10 4 for Eth = 100 eV. Smaller binding energies [upper line in Fig. [8b] result in up 
to two order of magnitude bigger rates but for threshold energies as low as -Eth 5 eV. 

Finally, in Table 2 we calculate the effects of time modulation and present the differential event 
rate for four different values of E' e in the case of massive mediator with mx = 1 GeV. We assume 
also a WIMP mass m x = 100 GeV and Z = 1. The H ratio is constant around 10% independent of 
the energy and the WIMP mass. So the modulation h of the total rate is also going to be around 
10 — 17%, which means that the difference between the maximum (here always in June 3nd) and 
the minimum (here always in December) is 20 — 34%. 

4.4 Experiment : The prospects of detecting single ultra low energy electrons 

As discussed in a previous section observation of light X- boson would require detectors with sub- 
keV sensitivities. The development of such detectors, having a low energy threshold and low noise, 
remains generally a daunting challenge for present-day and future low- background experiments. 

As shown in Fig. 6 the signal of low energy electrons produced by elastic collision process 
exhibits a maximum at energies around or even lower than 10 eV. At such energies a detector with 
single electron sensitivity will be required to reach a reasonable efficiency. The last ten years a 
particular effort is going on to develop ultra low threshold detectors in order to address low energy 
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neutrino physics [43-46]. This has been achieved for low mass detectors. We are, however, seeking 
an even lower energy threshold. 

Usual solid state detectors employed for dark matter projects have typical thresholds of a few 
keV. It is very difficult to combine sub-keV and big mass at the same time. For instance Ultra- 
Low-Energy Germanium detectors [36] are able to reach a threshold of a few hundred eV's, but 
they are limited to a modular mass of a few grams. Anyway the achieved energy threshold is still 
below our requirements. 

Single electron efficiency is achieved using detectors reaching very-high gains in order to cope 
with electronic noise. Gaseous detectors are good candidates. In such detectors high gains may eas- 
ily be achieved. Having been conceived as a TPC Micromegas detector (MS) [47], it is compatible 
with large drift volumes and operation at high pressure, an example of which are the HELLAZ [48] 
prototypes. A great advantage of this detector is the versatility of target material: various gases 
from the lightest (H2) to heaviest (Xe) could be used offering a large choice. 

One idea to increase the mass of the target material is to use the recently developed Spherical 
Proportional Counter (SPC). This detector consists of large spherical gas volume with central 
electrode and radial electric field. Charges deposited in the drift volume are drifting to the central 
sensor where are amplified and collected. A novel concept of a proportional sensor, a metallic 
ball having a radius of about 15 mm, located at the center of curvature, acting as a proportional 
amplification structure is used. It allows to reach high gas gains (> 104) and operates from low to 
high gas pressure. At such gains, provided the low electronic noise of this detector, single electron 
efficiency is easily achieved. 

The main advantages of the new structure relevant to our project are: 

• Simplicity of the design. 

• A single channel is used to read-out a large volume. 

• Robustness 

• The depth of the interaction, related to the rise time of the signal, is measured. This is 
important to apply fiducial cuts for background rejection purpose. 

• Low detector capacity < 0.1 pF, independent of the vessel size, allows very-low electronic 
noise, which is a key point toward achieving low energy threshold. 

• Versatility of the target material and density; the detector is compatible with a large variety 
of gases and could operate from low pressure to high pressure. This could be a precious tool 
to identify a possible signal out of backgrounds. 

A main concern of the proposed detection scheme is the minimal background level that will 
be reached by our system. By this one means that detector body and appropriate shield will 
be built with materials which are screened for low levels of natural and man-made radioactive 
impurities. Ordinary construction and shielding materials, however, do contain trace amounts of 
naturally occurring and man-made radionuclides which result in elevated background level; we need 
to design and fabricate the detector by careful material selection made out of low level activity. 

Unfortunately, however, there exists very little experience at the very low energy (sub keV) 
region where our detector will be operating. An example is a low background gaseous detector 
with sub KeV energy threshold developed for solar axion search [50] ; the reached background level 
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is quite low and is flat in the sub KeV energy range down to 250 eV. Our purpose is to further 
decrease the energy threshold down to about 10 eV. This region has never been explored and 
therefore reaching the desired low level activity becomes a new experimental challenge. Single 
electron backgrounds could be emitted by materials pulled by the electric field through thermionic 
emission. The advantage of the spherical detector is that at the external vessel the electric field is 
extremely low and therefore highly reduced thermionic emission is expected. 

The present prototype having a volume of lm 3 , filled with a gas at high pressure with a target 
mass of the order of 10 kg could fulfills sensitivity requirements for our project. We will search 
appropriate molecular gases having low binding energies and compatible with operation in the 
Spherical Proportional Counter detector [51]. 

At present it looks realistic to soon have a sphere of radius of 5 meters, which can be under 
a pressure of 5 bars. Thus, if one fills it with 80% Ar and 20% Isobutane (C±H\q), one can have 
212 Kg of Hydrogen. With this much Hydrogen using eq. (|4.68|) and a threshold of ~10 eV, we 
expect around 200 events per year for the parameters in (|4.67p . In models [19] where the mediator 
mass is very low, e.g. mx ~ 1 MeV, we expect an increase of the event rate by almost six orders 
of magnitude. Therefore, if a low energy experiment will be built it would possibly set the best 
limits on these kind of models. 

5 Conclusions 

Recent cosmic ray results from PAMELA, HESS and FERMI collaborations show an unexpected 
rising of positron events with energy that may be due to Dark Matter particle annihilations in the 
halo of our Galaxy. This Dark Matter particle "sees" the SM ones only through its interactions 
with an X-boson that couples to the SM gauge sector. Depending on the model, the mediator can 
be massless or massive with different couplings. In this article we study direct detection of this 
secluded type of dark matter employing nucleons or electrons with main emphasis in the latter 
case. 

Due to the small momentum transfei0 the massless case results in a large number of events 
that should have been seen by current nucleon recoiling direct detection experiments and therefore 
strong bounds on mixing parameters and couplings exist. Our work emphasizes the role of the low 
energy electron recoil in direct detection experiments and proposes a novel experimental avenue 
on how to proceed in searching for such low energy electrons. For simple hydrogenic atoms, and 
at low energy, E' e 10 eV, the cross section is enhanced by order of magnitudes compared to 
KeV recoil energies. In the neighborhood of low energies, the results depend highly on the binding 
energy of the ejected electron: the more loose the electron is the bigger the event rate becomes as 
expected. In this regard we considered two possibilities: 

• The process is mediated by the massive mediator X (our model III). 

In this case we do not have scattering off hadrons at tree level. So we do not have dominant 
constraints on the parameters of the model coming from the ongoing WIMP searches. Using 
the parameters of eq. (14.670 we have obtained fairly large cross sections for a Dirac WIMP. 
Employing the spherical TPC detector described above with a radius of 5 m under pressure 
of 5 Atm we have found that we could have about 200 counts in a year, assuming a threshold 
of 10 eV. It is possible, however, that our choice of parameters is a bit optimistic and we may 

7 For nucleons, the momentum transfer is ~ 2 MeV and energy transfer is ~ 2 KeV, while for low energy electron 
recoils they are ~ 50 KeV and ~ 10 eV, respectively. 
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have not considered all available constraints. Our results are also applicable to model-I. In 
this case however, due to the fact that couplings of the X-boson to hadrons appear at tree 
level, there exist strong constraints on the mixing parameter already from the nucleon direct 
searches [see eq. (|3.38j) ] . 

• The process is accommodated by the massless mediator (leptophylic version of model II) 
This mechanism is similar to that involving hadrons in section 3, one simply replaces the 
quarks by leptons. In this case we have found that the most stringent constraints on the 
parameters come from the standard WIMP searches. Thus using the parameters of Eq. 
(4.65) we have obtained with the above detector hundreds of events per year even with a 
(reduction) mixing coupling constant as low as k = 10- 10 for a Dirac WIMP. Such a huge 
signal cannot be seen by current experiments either due to lack of low energy threshold or 
because, experiments, like CDMS and XENON, are keeping only nuclear recoil events. We 
were surprised to find so large cross section. We now understand it, however, to be due to 
the photon propagator (1/q 2 ) 2 , which is favored by the fact that the momentum transfer 
is very low in the case of electrons. We should mention that, since the initial electron is 
bound, there is no infrared divergence and no need for a low energy cut off. It should be also 
noted that quark couplings to X-boson will come back through loop corrections even if they 
are forbidden at tree level by a symmetry which is eventually broken [see footnote 2]. Then 
current nucleon recoil experiments will be as important [see eq. (|3.32p ] and complementary 
to the electron ones. 

The above conclusions assume that the WIMP is a Dirac particle. If the WIMP is Majorana 
particle the rates are suppressed by approximately a factor (3 2 rj 10~ 6 . For both the above cases, 
annual time modulation effects are of the order of 20-30%, important enough to be noticed. 

We have limited the discussion of the rates in the case of hydrogen, since our cross section was 
evaluated using hydrogenic wave functions. Certainly the obtained rates will increase, if one can 
exploit the other atomic electrons with smaller binding energy. This situation was made manifest 
in our work with a judicious change of the binding energy [see Figs. 6b, 7c, 8b]. But then one should 
employ realistic wave functions. 

In a similar fashion one can treat other dark matter candidates like right handed neutrinos, 
which arise in models in which the ordinary Dirac type mass is forbidden due to a discreet symmetry, 
but communication with the leptons is allowed via exotic scalars [52-54] with masses in the 50 GeV 
region. It may also apply to other models involving exotic fermions and scalars recently proposed 
and reviewed in ref. [55]. 
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